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We show that the resolution "per absorbed particle" of standard absorption tomography can be 
outperformed by a simple interferometric setup, provided that the different levels of "gray" in the 
sample are not uniformly distributed. The technique hinges upon the quantum Zeno effect and has 
been tested in numerical simulations. The scheme we propose could be implemented in experiments 
with UV-light, neutrons or X-rays. 
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I. INTRODUCTION 

Absorption tomography is an important experimen- 
tal technique revealing the internal structure of mate- 
rial bodies. By measuring the attenuation of a beam 
of particles passing through a sample one infers the ab- 
sorption coefficient (density) of the sample in the beam 
section. The possibility of distinguishing two slightly dif- 
ferent densities of the material is often of vital impor- 
tance. Under otherwise ideal conditions the shot noise 
associated with the discrete character of the illuminating 
beam sets an upper limit to the resolution of absorption 
tomography: for instance, the shadow cast by a brain 
tumor might become totally lost in the noisy data. One 
possibility to overcome the fluctuations is to increase the 
intensity of the beam. However, in many situations, like 
in medicine for example, the intensity of the illuminating 
beam cannot be made arbitrarily high due to the damage 
provoked by the absorbed radiation. 

A significant step towards an "absorption-free tomog- 
raphy" came from quantum theory. It was demonstrated, 
both theoretically Q, 0] and experimentally , that to- 
tally transmitting and absorbing bodies can be distin- 
guished without absorbing any particles, by using an 
interferometric setup. This idea is in fact a clever im- 
plementation of the quantum Zeno effect jij and hinges 
upon the notion of "interaction- free" measurement [5(. 
A classical measuring apparatus (here the black sample) , 
placed in one arm of the interferometer, projects the il- 
luminating particle into the other arm, destroying inter- 
ference, freezing the evolution and forcing the particle to 
exit through a different channel from that it would have 
chosen had both arms been transparent (white sample). 

In practical applications, however, samples are nor- 
mally neither black nor white: they are gray. In this 
paper we endeavor to understand whether application of 
the quantum Zeno effect, which turns out to be ideal 
for discriminating black and white, might be advanta- 
geous also for the more practical task of discriminating 
two gray bodies with different transmission coefficients. 
More specifically, we ask: is it possible by quantum Zeno 
effect to reduce the number of absorbed particles while 



preserving the resolution? We show that this is indeed 
possible. Closely related questions have been recently in- 
vestigated by other authors 0, • Our conclusions are 
somewhat optimistic: we show that standard absorption 
tomography can be outperformed by a Zeno setup, pro- 
vided that the frequency of occurrences of the different 
levels of "gray" in the sample is not uniform. More to 
this, the Zeno setup, unlike the standard one, is endowed 
with two detection channels: as we shall see, this feature, 
if properly exploited, leads to even better performances 
in the Zeno case. 



II. QUANTUM ZENO EFFECT IN A 
MACH-ZEHNDER INTERFEROMETER 

We introduce notation and sketch the fundamental fea- 
tures of the quantum Zeno effect. We consider the Mach- 
Zehnder interferometric (MZI) scheme with feedback dis- 
played in Fig.^. A semitransparent object, whose trans- 
mission amplitude is r (assumed real for simplicity) is 
placed in the lower arm of the interferometer. The par- 
ticle is initially injected from the left, crosses the inter- 
ferometer L times and is finally detected by one of two 
detectors. The two semitransparent mirrors M are iden- 
tical and their amplitude transmission and reflection co- 
efficients are 



cos#l, s = siuOl (9l=k/4L), 



(1) 



respectively. Notice that both coefficients depend on L, 
the number of "loops" in the MZI. 

The incoming state of the particle (coming from the 
source at initial time) is 



|in) = 



(2) 



and we call "Zeno" and "orthogonal" channels the ex- 
traordinary (J) and ordinary (°) channels of the MZI, 
respectively. The total effect of the interferometer is 



V T = BA T B, B = 



A T = 
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FIG. 1: a) Scheme of the Zeno interferometric setup, b) 
Standard transmission experiment. S - source; M - semi- 
transparent mirror; o - orthogonal channel; z -Zeno channel; 
D - detector. 



In general, 

B = exp(-i9 L a 2 ), BB^ = B ] B = l, 



(4) 



where o~2 is the second Pauli matrix, while A T is not 
unitary (if r < 1 there is a probability loss). The final 
state, after the particle has gone through L loops, reads 



|out) = y T L |in) = (BA T B) L \in). 



A. White sample 



(5) 



The choice of the angle 9l in JU is motivated by our 
requirement that if r = 1 ("white" sample, i.e. no semi- 
transparent object in the MZI) the particle ends up in 
the "orthogonal" channel: 



VX, = B 



r>2L _ —i2L9 L a 2 _ —iira 2 /2 _ 



so that 



it)=V T L =1 \m} = 



(6) 



(7) 



This is easy to understand: each loop "rotates" the par- 
ticle's state by 29 l = ir/2L and after L loops the final 
state is "orthogonal" to the initial one (J2J. 



B. Black sample 

Let us now look at the case r = 0, corresponding to 
a completely opaque ("black") object in the MZI. We 
obtain 



V L 

V T=0 



B{A B 2 ) L B- 1 
B cos L 29 L 



1 -tan20 L 




L — ►oc 



1 






V r=0 



B- 1 



This yields QZE: 

| out) = V T=0 |in) = in) = 



(8) 



(9) 



In the infinite L limit the initial state is "frozen" and the 
particle ends up in the Zeno channel. 



C. Gray sample 



What happens if < r < 1? We easily get 



(l + r)c 2 



-sc(l + r) 



sc(l + t) t - (1 + t)s 2 



(10) 



The computation of V~ is straightforward but lengthy 
and yields a final expression which is elementary but 
complicated. However, we are mainly interested in the 
large-L limit, that for 



t < rf = (1 - sin7r/2L)/(l + sinvr/2L) (11) 



reads 



1 _ IlT^7 °( L r ) ) - () ( [ ^ , (|-)) 
0(L- X ) r^l + OiL- 1 )} 1+1 h 1 ' 



This is an interesting result: indeed 



< T < 1 



(13) 



analogously to JHJ . This shows that even for a semitrans- 
parent object, with transmission coefficient t^I. a bona 
fide QZE takes place and the particle ends up in the Zeno 
channel with probability one: 



|out) = V T |in) = |in), 



(14) 



III. 



DISTINGUISHING DIFFERENT SHADES 
OF GRAY 



A question arises : is it possible to distinguish differ- 
ent values of r (different "shades" or "levels" of gray) by 
the technique outlined above? This is not a simple task, 
for after a large number of loops L the particle ends up 
in the orthogonal channel only if r = 1 [see Eq. Q)]; by 
contrast, for any value of r ^ 1, the particle ends up in 
the Zeno channel [see (|14|l ] irrespectively of the partic- 
ular value of r. However, the asymptotic correction in 
the (1,1) element of in (|12fl is r dependent: the de- 
tails of the convergence to the limit l|13l) - (|14fl depend on 
the grayness of the sample. By exploiting this feature, 
we shall now show that it is indeed possible to resolve 
different gray levels by QZE, within a given statistical 
accuracy. 

We start by observing that if one performs a standard 
transmission experiment, by shining a particle beam on 
a semitransparent object in order to measure the trans- 
mission coefficient r, see Fig. the detection and ab- 
sorption probabilities read 



P' d (r) 



Pair) 



(15) 
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The statistics is binomial. 

On the other hand, if one uses the Zeno configuration 
sketched in Fig.^,, the final state of the particle after L 
loops in the MZI is, from (JT35J, 



u z 

U n 



7T 2 1 + T 

8L 1-r " 

0(L" 



0(L- 



(16) 



where u z and u Q are the amplitudes in the Zeno and 
orthogonal channels, respectively. Both these quantities 
are real. Let (0 < t < 1) 



p z (r) = ^ = i-ZLl±l + (L- 2 ), 



Po(t) 



4L 1 - r 
0(L- 2 ), 



(17) 



p a (r) = l- Pz(T )- Po (r) = | E i±^ + 0(L- 2 ) 

be the probabilities that the particle is detected in the 
Zeno, orthogonal channel or is absorbed by the semi- 
transparent object, respectively. We assume that a fixed 
number of particles N is sent in the MZI during an exper- 
imental run. In this situation the distribution of particles 
in the Zeno, orthogonal or absorption channels follows a 
trinomial statistics with probabilities 1171) . 

It seems natural to think that, since by increasing the 
number of loops L, p a vanishes much faster than p a , for 
large L the distribution is practically binomial with p a + 
Pz ~ 1 However, as we shall see, the presence of 

a small trinomial component, p , will play an important 
role, enabling the Zeno method to perform much better 
than the standard one. 



IV. AN INTRINSIC LIMIT FOR BINOMIAL 
STATISTICS: THE CRAMER-RAO BOUND 

We are now ready to discuss the possibility of a "Zeno 
tomography." The goal is to get information about the 
distribution of the absorption coefficient in the sample, 
absorbing as few particles as possible. We will accom- 
plish this in two steps. First, using estimation theory, we 
show that if one limits one 's attention to binomial statis- 
tics, the Zeno estimation of any level of gray of one pixel 
(i.e. r continuously distributed between and 1) cannot 
perform better than the standard method. At best, both 
methods are equivalent. This is bad news. However, the 
very proof of the above-mentioned statement will show 
that there are two ways out: first, we will see that Zeno 
performs better when one wants to distinguish two levels 
of gray that are not equally populated in the sample (this 
requires some prior knowledge about the distribution of 
grays in the sample). This is good news, for it enables 
one to find a method that in some cases works better 
than the standard one. Second, one is led to think that 
the introduction and exploitation of a trinomial statistics 
can enable the Zeno method to perform better. 

Let us start from the estimation of any level of gray. 
In this case one tries to estimate r 2 from the counted 



number of particles, absorbing as few particles as possible 
for the requested precision. To perform this task in an 
optimal way one should find optimal estimators for each 
scheme. A lower bound on the variance of an unbiased 
estimator T of the parameter T (here T = t 2 ) is the 
Cramer-Rao lower bound (CRLB) 0, 



(AT) 2 >- = 



d 

— lnp(n|T) 



1 2\ 



(18) 



where F is the Fisher information, p(n\T) the probabil- 
ity of observing n particles conditioned by the value T of 
the unknown parameter and (...) denotes ensemble av- 
erage with respect to n. The probability p is binomially 
distributed in both the standard case and the Zeno case, 
which yields 



(AT st ) 2 
(Afzc) 2 



> 



> 



r 2 (l 



N 

4t 2 (1-t) 3 (1 + t)L 



(19) 



for standard and Zeno tomography, respectively, N being 
the (fixed) number of input particles in both cases. Ex- 
pressing the above inequalities l|19l) in terms of the num- 
ber of absorbed particles — Np' a and — Np a , 
they both reduce to the same bound 



AT opt > 



-(1 -r 2 ) 

jy(st),(Ze) 



(20) 



showing that the CRLB's for standard and Zeno tomog- 
raphy are the same, given the number of absorbed par- 
ticles. Further, it is trivial to show that the unbiased 
estimator given by the relative frequency of transmit- 
ted particles, T st = n t /N, saturates the CRLB l)20fl. 
Hence, if one neglects the output of the ordinary chan- 
nel p and considers the statistics (I17f) practically bino- 
mial, the Zeno estimation can be at most as good as the 
standard one: it cannot be better. 



V. TWO STATISTICAL PROTOCOLS 

A. Binomial (single-channel) protocol 

In spite of the conclusions of Sec. lIVI we will now con- 
struct a protocol and show that by QZE one can achieve 
a resolution that is superior to the "ordinary" resolution 
obtained in a standard transmission experiment. Notice 
that p Q (r) in (|17l) . unlike p' a (r) in i|15|) . is an increasing 
function of r. Therefore, with respect to absorbed par- 
ticles, the Zeno tomographic image (for sufficiently large 
L) yields a kind of negative of the standard absorption 
tomographic image. This can be given a rather intuitive 
explanation: indeed, the absorption probability in \TJ$ 
reduces to the same form as the standard one (|15[) . i.e. 



Pa(r) 



(21) 
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by introducing an "effective" transmission coefficient 

T-eff = VT^pI. (22) 

For example, if we take t% = 0.98, t 2 = 0.99 and 
choose L = 12000, then, according to Eq. I|17|l . we get 
w 0.99 and r c z ^ 2 w 0.98. The two gray levels are 
interchanged by the Zeno apparatus. If most of the sam- 
ple has transmission coefficient t 2 the absorbed energy is 
reduced by using the Zeno setup. 

A more precise comparison of the performances of the 
Zeno and standard techniques can be given in the frame- 
work of decision theory. For simplicity let us focus on 
distinguishing only two gray levels t\ and t 2 {t\ < t 2 ) 
corresponding to hypotheses H\ and H 2 that occur in 
the sample with frequencies 

P (H 1 ) = a, P (H 2 ) = l-a. (23) 

With this simplification we lose no generality since the 
tomography with M gray levels can always be split into a 
sequence of pairwise decisions between two adjacent gray 
levels. 

We will proceed in two steps. First we will assume that 
in the Zeno configuration of Fig. ^ all output particles 
are collected at a single detector. In other words, the 
Zeno and orthogonal channels are considered as a single 
output and the statistics l|17l) is binomial (p z +p a ,p a ). 
Each particle is then either absorbed or transmitted (and 
detected) by the Zeno apparatus. Obviously, by merg- 
ing the two output channels together, some information 
about the sample is wasted. We know, however, that 
this strategy will be optimal if the number of loops L 
is very large. Then there are almost no particles exit- 
ing via the ordinary channel n a « which can then be 
safely ignored. A better and more general strategy will 
be studied later. 

Since both experiments obey the same (binomial) 
statistics we use the notation of the Zeno experiment. 
The analysis of the standard experiment is similar. If no 
distinction is made between the Zeno and ordinary chan- 
nels, the decision is based on the number of absorbed 
particles n a . If n a is smaller than or equal to a decision 
level n„, then Hi is chosen; otherwise H 2 is chosen. The 
probability of making an error in identifying the gray 
level of a given pixel is 

P e = aP(H 2 \H x ) + (l-a)P(H x \H 2 ), (24) 

where 

P{H 1 \H 2 )= P(n a \H 2 ) (25) 

is the probability of choosing Hi when Hi is true and 

p(n a \H 2 ) = ^p a (T 2 r«{l-p a (T 2 )] N - n ° (26) 

is the binomial probability of absorbing n a particles when 
H 2 is true. [For P(H 2 \H\) the summation is over n a > 



n d ] . An optimal protocol is given by determining n d a that 
minimizes the error (|24f) . 

Alternatively, one defines the likelihood ratio 

£(Ti|iV,rt a ,a) 

H= T? — iTt i V ( 27 J 

£{t 2 \N, n a ,l - a) 

where 

C{n\N,n a ,a) = Pb(Hi)p(n a \Hi) 

= a (*)p a (T i )'*[l-p a (Tj] N - n ' (28) 

represents the likelihood of hypothesis Hi (i = 1,2). The 
optimum decision level n d is determined by solving for 
equal likelihoods 

R = l. (29) 

In both cases one gets 



log 


[^]- 


- ATlog 


T- 
l- 


-Pa(Tl) 
-Pa(T 2 )_ 


log 


Pa(r 2 )_ 


- log 


"l- 

l- 


Pa(ri) 

Pa(T 2 ) _ 



(30) 



and, substituting in Il24(l . 

P e = a{l - Bj[N- nil + nt;l- pain)]} 

+ (1 - a)Bj[N - h d a , 1 + r~i d a ; 1 - p a (r 2 )}, (31) 

where Dj(a, b: z) is the regularized incomplete Beta func- 
tion 01 an d h d is the greatest integer less than or equal 
to n d a . The mean number of absorbed particles is 

JVf*) = N[ap a ( Tl ) + (1 - a)p a (r 2 )}. (32) 

By plugging Eqs. 1(201) and O in Eq. (JSTJ, the average 
probability of error l|31|l can be expressed as a function 
of a, t = t\, dr = t 2 — T\ and ■ The probability of 
error for the standard setup is obtained in a completely 
analogous way. 

The performances of the Zeno and standard methods 
are compared in Fig. |5J First, the (exact) Lth power of 
the matrix V T in IjlOfl is evaluated, then for a given error 
rate P e , the number of absorbed particles is calculated 
by solving numerically Eq. (|31|l . Their ratio is shown as 
a function of a for a few values of the transmission co- 
efficient r. Notice that the exposition of the sample can 
be significantly reduced if the distribution of gray lev- 
els in the sample is not uniform. For instance, a reduc- 
tion factor of 2.5 is obtained when the sample consists 
of 97% of dense material and 3% of the less absorbing 
one, a — 0.97. Such parameters are typical for struc- 
tural analyses: a small structural defect (crack) inside a 
thin sample would typically show small contrast (St « 1) 
with the surrounding almost transparent (t ~ 1) mate- 
rial, while its area would be small compared to the area 
of the sample (a w 1). 
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FIG. 2: Ratio of the number of absorbed particles in the 



Zeno (Na ) and standard (AT„ ; ) setup. The smaller the 
ratio in the graph, the less irradiation in the Zeno appa- 
ratus (for the same resolution). P e = 0.5%; dr — 0.02; 
r = {0.8, 0.9, 0.95, 0.97}; L = 2000. 



B. Trinomial (two-channel) protocol 

The binomial decision strategy outlined in the previous 
subsection is not the optimum one. Unavoidable losses 
and other imperfections of real experimental devices set 
a strict limit on the maximum number of loops that can 
be achieved in a laboratory. In such case the ordinary 
channel can no longer be ignored. The data consist then 
of the two component vector (n z ,n ) of the numbers of 
particles counted in the Zeno and ordinary output chan- 
nels. The decision will be based on both these numbers. 

The decision levels are readily obtained from the equal 
likelihood criterion 



R 



£{Ti\N,n z ,n ,a) 
C(T 2 \N,n z ,n , 1- a) 



1, 



(33) 



where 



C(t\N, n z , n Q , a) 



aN\ 



n z \n \n a \ 



p z (TT> Po (TT°p a (Ty 



This equation is to be solved for the decision vector 
(n d ,n d ). Equation ffify is just one condition for the two 
unknowns n z and n Q , so there exists a one-parametric 
family of solutions. By plugging (|31Jl into one easily 
finds 



n d z - a{T 1 ,T 2 )n d = 6(ti,t 2 , a), 
where the coefficients a and b read 
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(35) 
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FIG. 3: Typical decision levels of the binomial (dashed lines) 
and trinomial (solid lines) decision strategies when awl. 
This figure corresponds to the simulation shown in the last 
row of Fig. [KJi, where M = 3, the triangle is divided into 
M — 3 regions and the three gray shades are labelled as white, 
gray and black, respectively. 



In a two-dimensional representation each possible exper- 
imental outcome (n 2 ,n ) is represented by a point lying 
inside the triangle {0 < n + n z < N} shown in Fig. [3] 
Equation (|35J) divides this triangle into two regions. All 
experimental outcomes that fall within the same region 
issue the same decision. In the general case of M different 
gray levels there are M — 1 equations (|35() defining M — 1 
in general non-parallel lines dividing the square into M 
strip-like regions. This is shown in Fig. for M = 3. 

An interesting situation arises when the coefficient a 
in H35|) becomes close to unity. In that case, the decision 
level is the line n z — n d Q = const. Let us recall that the 
decision levels of the binomial decision strategy discussed 
in the previous subsection were n d = const, (see l|5U)l). or, 
equivalently, n z + n d Q = const. Hence if L, T\ and t 2 are 
such that a ~ 1, the decision levels of the binomial and 
trinomial decision strategies are orthogonal to each other. 
This is shown in Fig. |3J Under such conditions one can 
expect further gain in the precision of the Zeno apparatus 
as compared to standard absorption tomography. This 
regime was chosen for our computer simulations of the 
following section. 

Notice that the steepness of the decision lines (|35fl de- 
pends only on the absorption of the corresponding adja- 
cent gray levels. It depends neither on their frequencies, 
nor on the total number of incident particles. 

Finally, let us discuss the limit L — > oo. When the 
number of loops L increases, p — > much faster than 
p a [see Q170 ] and one is allowed to put n d = and n d = 
N - 



in Eq. I|35|) . which reduces to 
N - 6(ri,r 2 ,a). 



(38) 
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FIG. 4: The object to be reconstructed: a cell of Giardia 
lamblia, one of the most primitive eukaryotes. The original 
picture has been reduced for simplicity to three levels of gray: 
white, gray and black, occurring with frequencies a m = 0.02, 
Oi g = 0.07 and Of, = 0.93 respectively. 

By substituting p z (r) = 1 — _p a ( r ) in Ip ffi -ffB jl . one reob- 
tains the binomial condition (|30|) : the binomial strategy 
becomes optimal in this limit. 



VI. SIMULATIONS 

We have seen that the Zeno technique can reduce the 
level of absorption without losing resolution (compared 
to the standard technique). (Alternatively, the Zeno 
setup can yield an improved resolution, while keeping the 
absorption at the same level of the standard setup.) The 
object in Fig. 01 is a cell of giardia lamblia, a protist, one 
of the most primitive eukaryotes. Giardia has been called 
a "missing link" in the evolution of eukaryotic cells from 
prokaryotic cells. The number of gray levels in the figure 
has been reduced to three to make the analysis simpler: 
white, gray and black, t w = 0.99, r g = 0.96, u = 0.8, 
occurring with frequencies a w — 0.02, a g — 0.07, and 
as = 0.93 respectively. Figure shows the results of a 
numerical simulation, performed with the standard and 
Zeno methods, the latter for L = 10 and L = 165, for dif- 
ferent numbers of absorbed particles N a . In each frame 
the standard and the two Zeno reconstructions are com- 
pared, together with the pixels that have been misinter- 
preted. Figure El confirms the expectation based on the 
asymptotic formulas Ijl7|) : in general, provided that the 
object contains a small fraction of more transparent pix- 
els and a larger fraction of more absorbing material, the 
Zeno setup yields a better resolution for a given irradia- 
tion. Clearly, a significant improvement with respect to 
standard absorption tomography is achieved for as few as 
L = 10 loops. The improvement is very large for L = 165. 

The number of absorbed particles increases from (a) 
to (d) in Fig. [SJ Observe that in (a) the standard re- 
construction fails completely, while the outline and basic 
shape of the object can be recognized already in the Zeno 
reconstruction with L = 10. In (c) the Zeno reconstruc- 
tions are quite good, while standard tomography does 
not detect white pixels in the object. When the intensity 



of the illuminating beam is increased further, in frame 
(d), all the reconstructed images become visually hard 
to tell from the sample, but the error rates of the Zeno 
apparatuses are still much better (by a factor three or 
more), as shown by the number of misinterpreted pixels. 

It is worth commenting on the distribution of misin- 
terpreted pixels. Clearly, in all the cases analyzed, it is 
not uniform. In general, when the distribution of gray 
levels in the sample is not uniform, any reconstruction 
technique tends to perform better in the "background," 
while making more mistakes in the region where the 
"structure" is present. The improvement due to the 
Zeno method becomes apparent if one looks in partic- 
ular at FiguresOJb) and (c): in these cases, interestingly, 
the standard method yields more mistakes in the back- 
ground; this is an unpleasant feature, if one is interested 
in detecting small irregular structures in a more or less 
uniform background. The features of the distribution of 
misinterpreted pixels require more careful study and their 
comprehension might lead to additional ideas. 

Any increase in the number of loops L in the interfer- 
ometer makes the difference between standard and Zeno 
tomography even bigger. Clearly, this is more demanding 
in terms of experimental realization. 



VII. CONCLUSIONS 

We have shown that a quantum Zeno tomography is 
possible and performs better than standard tomography 
if a given prior knowledge about the distribution of grays 
in the sample is available. This is a common situation in 
radiography, where one is often interested in detecting a 
small structure in a uniform background, like for instance 
in the analysis of small structural defects. 

In our numerical simulations we have illustrated some 
situations in which the resolution is improved by the Zeno 
method, for a given number of absorbed particles. Alter- 
natively, for a given resolution, the Zeno method per- 
forms better, absorbing less particles. This can be inter- 
esting in applications, for instance if one wants to limit 
the damage provoked by the absorption of radiation with- 
out losing in resolution. 

It is obvious from Fig. |2 that an even larger improve- 
ment is possible for almost transparent samples, provided 
that a is close to unity. This means that there is no fun- 
damental limit on the improvement that can be achieved 
over the standard setup: in other words, there is no "op- 
timal" configuration. 

There are additional issues that deserve careful study. 
For instance, the effects due to a Poissonian beam (total 
number of incoming particles N not fixed) and a complex 
transmission coefficient r. 

Let us also comment on experimental feasibility. Fig- 
ure [5] shows that an experimental test of the Zeno tomo- 
graphic technique should not be as difficult as one might 
think: simulations have been performed for as few as 
L = 10 loops in the interferometer, giving better results 
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FIG. 5: Comparison of standard and Zeno tomographic techniques. In each frame: top left=reconstruction by standard 
technique; top right=misinterpreted pixels by the standard technique; center left=reconstruction by Zeno technique with 
L = 10; center right=misinterpreted pixels by the Zeno technique with L = 10; bottom left=reconstruction by Zeno technique 
with L — 165; bottom right=misinterpreted pixels by the Zeno technique with L = 165. The mean number of absorbed particles 
per pixel (irradiation) is N a = 1.7, 2.3, 4 and 13 for frames (a), (b), (c) and (d) respectively. The total number of particles N 
(total energy) scales approximately like 3-/V a , 1.8N a , and 6.5iV a for top, center and bottom reconstructions, respectively. We 
used t w = 0.99, r fl = 0.96, and Tt = 0.8. The sample consists of 10,000 (= 100 x 100) pixels, where white, gray and black occur 
with frequencies a w = 0.02, a g — 0.07 and ai = 0.93, respectively. The number of misinterpreted pixels are (top to bottom): 
(a) 968, 786, 315; (b) 942, 596, 212; (c) 717, 382, 68; (d) 205, 69, 0. 



than the standard method. It is reasonable to think that 
a Zeno setup with a much larger number of loops can be 
built for UV light (highly absorbed by some biological 
samples). Also, by changing the light wavelength, one 
could efficiently "observe" different regions of the sample 
(or slightly different samples). Moreover, the experimen- 
tal configuration we have proposed (photons in a MZI, 
like in Fig. ^) is certainly not the only conceivable one. 
Phase imaging and tomograph y ha ve been demonstrated 
for both X-rays and neutrons [llj. More to thiSj Rauch 
and collaborators, with the VESTA apparatus [12|, have 
been able to keep neutrons in a one-meter long perfect 
crystal storage system ("resonator") for a few seconds, 
so that the neutrons bounce back and forth between two 



mirrors several thousands times. This would lead us to 
the full asymptotic (L 3> 1) regime considered in Fig. [21 
and last row of Fig. \5\ where the Zeno method can per- 
form much better. There is hopefully more to come. 
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